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The  purpose  of  this  note  is  to  state  the  solution  to  the  inhomogeneous  matrix 


differential  equation 
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I. 


TERMINOLOGY 


Suppose  that  we  are  given  the  time  varying  n  x  n  matrices  A(t),  B(t),  U(t).  We 
shall  assume  that 

a.  the  elements  of  A(t)  and  B(t)  are  continuous  functions 
of  the  time  t 

b.  the  elements  of  U(t)  are  piecewise  continuous  functions 
of  t. 

We  shall  seek  the  solution  of  the  matrix  differential  equation 

^  X(t)  =  A(t)  X(t)  +  X(t)  B(t)  +  U(t)  (1) 

subject  to  the  initial  condition 

X(t  )  =  X  (2) 

o  o 

where  X(t)  is  an  n  x  n  matrix. 

II.  THE  HOMOGENEOUS  CASE 

Bellman  in  Reference  [  1]  (page  175)  considers  the  homogeneous  equation 

X(t)  =  A(t)  X(t)  +  X(t)  B(t)  (3) 

subject  to  the  initial  condition 

X(t  )  =  X  (4) 

o  o 

His  result  is  that  the  solution  of  (3)  is  given  by  the  relation 

X(t)  =  4>(t;to)  Xq  *(t;tQ)  (5) 

where  4>(t;to)  is  a  nonsingular  fundamental  matrix  which  satisfies  the  differential 
equation 
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(6) 


—  $(t;t  )  =  A(t)  $(t;t  );  4>(t  ;t  )  =  I 
at  o  o  o  o 

and  where  \t(t;t  )  is  a  nonsingular  fundamental  matrix  which  satisfies  the  differential 
o 

equation 


-T-  *(t;t  )  =  *(t;t  )  B(t);  *ft  ;t  )  =  I 
dt  o  o  o  o 


(7) 


III.  THE  INHOMOGENEOUS  CASE 

We  claim  that  the  solution  of  the  differential  equation 


^  X(t)  =  Aft)  Xft)  +  X(t)  Bft)  +  Uft) 


(8) 


with  Xft  )  =  X  is  given  by 


o 


Xft)  =  $ft;to) 


X  +  f  $  1  (r;t  )  U(t)  *-1  (r,t  )  dr 
o  J  t  o  o 


(9) 


To  see  this,  differentiate  (9)  with  respect  to  t  to  obtain  (we  use  dots  to  indicate 
differentiation)  the  relations 


Xft)  =  $(t;to) 


Xo  +  j  t  (T»t0)  U(t)  ^_1(r,to)dT 

o 


*<«o> 


+  *ft;t  )  1  (t;t  )  Uft)  )  *(Gt  ) 

o  o  o  o 


+  ^(t;tQ) 


X  +  f  (r,t  )  U(i)  *_1  (T;t  )dr 

o  J  t  o  o 


*ft;t  ) 

O 


=  Aft)  *ft;t  ) 

[x  +  r 

o 

[  O  J 

t 

0 

xft) 


-  $(t;t  )|  X  +  [  s'^t  )U(r)^_1  (r,t  )dr]  ^(t;t  )B(t)+U(t) 
olojt  o  o  0 

- - - - Q - - - i _ _ 


Xft) 


(10) 
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and,  so, 


X(t)  =  A(t)  X(t)  +  X(t)  B(t)  +  U(t) 


IV.  TIME  INVARIANT  CASE 


If  A  and  B  are  constant  matrices,  then 

A(t-t  ) 

4>(t;tQ)  =  e 


^(t;to)  =  e 


(H) 


(12) 


(13) 


and,  so,  the  solution  reduces  to 


X(t)  =  e 


pt  -A(T-t  )  -B(T-t  ) 

X  +  \  e  °  U(r)e  °  dT 
o  J  t 


B<'-to> 

e 


(14) 
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